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ABSTRACT 

We describe how AGN-produced cosmic rays form large X-ray cavities and radio lobes in the hot 
diffuse gas in galaxy groups and clusters. Cosmic rays are assumed to be produced in a small shocked 
region near the cavity center, such as at the working surface of a radio jet. The coupled equations 
for gasdynamics and cosmic ray diffusion are solved with various assumptions about the diffusion 
coefficient. To form large, long-lived cavities similar to those observed, the diffusion coefficient must 
not exceed k ~ I0 28 cm 2 s _1 in the hot gas, very similar to values required in models of cosmic ray 
diffusion in the Milky Way. When k 10 28 cm 2 s -1 , cosmic rays are confined within the cavities for 
times comparable to the cavity buoyancy time, as implied by observations of X-ray cavities and their 
radio synchrotron emission. Collisions of proton cosmic rays with thermal plasma nuclei followed by 
7T° decays can result in enhanced gamma ray emission from the cavity walls. 

Subject headings: X-rays: galaxies - galaxies: clusters: general - X-rays: galaxies: clusters - galaxies: 
cooling flows 
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1. INTRODUCTION 

In recent years high resolution X-ray observations of 
the hot gas in galaxies, groups and clusters have re- 
vealed cavities or bubbles having sizes of a few kilopar- 
secs. About 20 percent of all groups/clusters are cur- 
rently known to contain these cavities (e.g. Bohringer et 
al. 1993; McNamara et al. 2000; Birzan et al. 2004), but 
they may appear intermittently in all clusters. In previ- 
ous computational studies X-ray cavities have been pro- 
duced by the introduction of ultrahot gas that expands, 
displacing the ambient hot gas in a purely gasdynamical 
manner (e.g. Churazov et al. 2001, 2002; Quilis et al. 
2001; Briiggen & Kaiser 2002; Briiggen 2003; Robinson et 
al. 2004; Jones & De Young 2005). However, many cav- 
ities are filled with radio synchrotron and inverse Comp- 
ton emission produced by relativistic electrons. Although 
it is often claimed that cavities are inflated by the pres- 
sure of cosmic rays, cavity formation by cosmic rays has 
not been investigated in detail. The interaction of cos- 
mic rays with hot gas differs from ordinary gas dynamics 
because of the important role of cosmic ray diffusion. In 
the discussion below we describe the combined dynam- 
ics of gas and diffusing cosmic rays to determine those 
properties of cosmic rays - energy density, luminosity, 
diffusion coefficient, etc. - that are best suited to form 
X-ray cavities similar to those observed. 

It is widely assumed that mass accretion onto a cen- 
tral massive black hole (i.e. AGN) is the ultimate energy 
source that forms these cavities, but intermediate details 
are lacking. We assume here that cosmic rays are the 
agency that transmits AGN energy to the gas. Cosmic 
rays can be formed in shock waves at the tip (working 
surface) of a radio jet and this is supported by many ob- 
servations. Cosmic rays may also be produced by shocks 
or other processes near the central black hole. In cither 
case, as cosmic rays propagate into the hot ionized gas, 
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they are scattered by irregularities in the magnetic field 
which, if frozen into the hot gas, convey cosmic ray mo- 
mentum to the gas. 

Once a low-density X-ray cavity forms at some radius 
in the atmosphere of hot cluster gas, its buoyancy trans- 
ports it radially outward. Therefore, to be visible as iso- 
lated, coherent features, X-ray cavities must form during 
time scales t; f, e that are considerably less than the buoy- 
ancy time tbuoy, which is in turn somewhat longer than 
the dynamical (freefall) time in the local galactic or clus- 
ter potential. However, cavities must not be formed ex- 
plosively on very short time scales since this would create 
powerful shocks, strongly heating the gas in the cavity 
rims that is not observed. Indeed, the gas in cavity rims 
is often the coolest gas in the entire hot gaseous atmo- 
sphere (e.g. Fabian 2001, 2002; Brighenti & Mathews 
2002). 

We present below a series of spherical, one- 
dimensional, time-dependent calculations that describe 
how cosmic rays create X-ray cavities. For simplicity we 
assume that the gas has uniform density and tempera- 
ture prior to the introduction of cosmic rays. In this first 
study we concentrate on the formation and duration of 
cavities and ignore the additional complication of grav- 
ity which however is essential to understand the buoyant 
phase. 

2. EQUATIONS AND COMPUTATIONAL PROCEDURE 

We assume that the cosmic rays (CRs) are produced 
in a shock located in a volume that is small compared to 
the cavity radius. The CR-gas interaction can be most 
simply expressed with four Eulerian equations: 
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ft +V 
^ + (u-V)u) =-V(7>- 



Pc) 



de 
dt 



V • ue = -P(V • u) 



(2) 



(3) 



de 

— ^ + V-ue c = -P c (V-u) + V-( K Ve c ) (4) 



2 



where artificial viscosity terms are suppressed. A sim- 
ilar set of equations has been discussed by Jones & 
Kang (1990). Pressures and thermal energy densities 
in the plasma and cosmic rays are related respectively 
by P = (7 — l)e and P c = — l)e c where we as- 
sume 7 = 5/3 and -f c = 4/3. The cosmic ray dynam- 
ics are described by e c , the integrated energy density 
over the cosmic ray energy or momentum distribution, 
e c oc / EN(E)dE oc / p A f(p){l+p 2 )- 1 / 2 dp. In general 
the CR diffusion coefficient k may depend on the particle 
momentum, n{p) ocp 5 " 6 , as suggested by CR models 
in the Milky Way (Malkov 1999), but we do not consider 
this possibility here and regard both e c and k as mean 
values integrated over the CR spectrum. The first three 
equations above describe the dynamical evolution of the 
hot gas influenced by the combined CR and gas pressure. 
The last equation describes the simultaneous gasdynam- 
ical advection and diffusion of the CR energy density. A 
mass conservation equation for the CRs is unnecessary 
because of their negligible rest mass. 

Cosmic rays and hot plasma are coupled by mutual 
interactions with a small (and otherwise dynamically in- 
significant) magnetic field with topological irregularities 
that scatter and diffuse cosmic ray particles on small spa- 
tial scales. However, no magnetic terms need explicitly 
appear in the equations provided magnetic stresses are 
small compared to P and P c . This assumption is justified 
by the small magnetic fields observed in group/cluster 
gas, ~ 1 — 10/iG (Govoni & Feretti 2004), corresponding 



to energy densities ~ B z 
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generally much less than the thermal energy density in 
the hot gas. Another implicit assumption in Equation (4) 
is that the (Alfven) velocity of the magnetic scatterers is 
small relative to the hot gas (e.g. Drury & Falle 1986). 
This assumption is also reasonable for our purposes since 

the Alfven velocity v A = B/(4tt / 9) 1 / 2 = 2rie 1/2 B(pG) 
km s _1 is generally very much less than the sound or 
flow speeds for magnetic fields and plasma electron den- 
sities of interest here. 

Our neglect of magnetic terms in Equations (l)-(4) also 
applies to the radio lobes. Observations of inverse Comp- 
ton (IC) X-rays (i.e. upscattered cosmic microwave pho- 
tons) from many radio lobes in X-ray clusters permit esti- 
mates of both the cosmic ray and magnetic energy densi- 
ties (e.g. Longair 1994), independent of the cquipartition 
assumption (Burbidgc 1959) which may not be physi- 
cally justified. Croston et al. (2005) observed keV IC 
emission in 33 FRII sources in cluster cooling flows and 
found magnetic fields within 0.3-1.3 of equipartition as- 
suming pure electron (e ) cosmic rays. The magnetic 
fields B ~ 10^iG in the radio lobes are comparable to 
those estimated from Faraday depolarization through- 
out all the hot gas (e.g. Govoni and Feretti 2004) and 
are dynamically unimportant. Croston et al. argue that 
cosmic ray protons are not a dominant component since 
cquipartition can be achieved with electrons alone. How- 
ever, Croston et al. do not constrain the radio lobes to 
be in pressure equilibrium with the local hot gas. In 
their study of radio emission from FRI lobes, Dunn & 
Fabian (2004) assumed pressure equilibrium with the sur- 
rounding gas and found that the pressure of cosmic ray 
electrons exceeded equipartition values in all sources, i.e. 
the particle energy density exceeds that of the magnetic 



field. About half of the lobes they studied require a dom- 
inant population of cosmic ray protons, in contrast with 
the results of Croston et al. For our purposes here it is 
not necessary to specify the CR composition, either elec- 
trons or protons can dominate. But we do assume that 
the total CR energy density is not substantially reduced 
during the cavity evolution time by losses due to syn- 
chrotron emission or interactions with ambient photons 
and thermal particles. 

In the model we adopt in Equations (l)-(4), we assume 
that the random fields that scatter CRs dominate over 
fields that are coherent on the scale of the radio lobe, con- 
sequently the diffusion coefficient k can be represented as 
a scalar rather than a tensor. This seems plausible be- 
cause the shapes of X-ray cavities are to first order spher- 
ical, not highly elongated. The CR diffusion described 
by Equations (l)-(4) are appropriate for CR-scattering 
field irregularities (Alfven waves) that are coupled to the 
inertia of the plasma. It is less clear if the CRs diffuse 
in the same manner relative to the field inside the cav- 
ities/lobes where the plasma density may be very low, 
but we assume that they do. 

3. LAGRANGIAN CALCULATIONS 

In this initial study of the combined dynamics of hot 
gas and cosmic rays we solve the one-dimensional spher- 
ical version of the equations above. While solutions can 
be found in either Eulerian or Lagrangian computational 
grids, we discuss the Lagrangian solutions here because 
they are conceptually simpler, require fewer computa- 
tional adjustments, and show more dramatically how the 
gas is cavitated by the cosmic rays. The Lagrangian 
equations are 
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de c 
dr 



(5) 
(6) 
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The first four equations are the usual Lagrangian gas- 
dynamical equations where m is the mass coordinate 
and e = e/p is the specific thermal energy. Differenced 
versions of these equations are solved explicitly at each 
timestep for time-advanced quantities - u, r, p, £ and e c 
- in the same order as the equations above. Equations 
(8) and (9) are solved term by term using operator split- 
ting. The diffusion term in Equation (9) is solved with 
an implicit Crank- Nicolson difference scheme, stable (but 
not necessarily accurate) for any time step. Although 
the Lagrangian grid necessarily follows the gas as it is 
pushed away from the source of cosmic rays, the cos- 
mic rays continue to fill the central cavity left behind. 
Consequently, the cosmic ray diffusion term in Equa- 
tion (9) is solved on a stationary Eulerian grid within 
the cavity and on a time-varying Lagrangian grid in the 
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gas. Both Eulerian and Lagrangian grids are updated 
as necessary at each time step. The computational time 
step At = amin(Ai c , At a) is a fraction a < 1 of the 
minimum of the Courant time step At c = Ar/c s (with 
sound speed c s = [(7P + J C P C )/ p] 1 ^ 2 ) and the time step 
for stability of explicit differencing of the diffusion term, 
Atd = {Ar) 2 /An. In zones where the gas is strongly 
compressed, Ar becomes small, slowing the calculation. 
Therefore we remove zones by combining adjacent La- 
grangian zones (conserving mass and energy) whenever 
Ar drops below an r-dependent minimum value. 

Since the hot gas is strongly compressed in a large 
region ahead of the diffusing cosmic rays, the artificial 
viscosity Q needs some modification to avoid excess dis- 
sipation. In addition to the standard Richtmeyer-von 
Neumann definition of Q which is nonzero only in in 
compressing zones, we require that (1) Q = if the gas 
gradient 91ogP/91ogr < 0.3 across a zone and (2) Q is 
not allowed to exceed the local gas pressure P. 

For these calculations we regard the CR diffusion co- 
efficient k as an unknown parameter. In the absence of 
detailed information about the nature of the magnetic 
scattering, it is difficult to derive k from first principles. 
However, k can be estimated for the X-ray cavity prob- 
lem from the global properties of observed cavities. As 
CRs diffuse away from a small source, their pressure gra- 
dient displaces the ambient hot plasma until a cavity of 
radius ri b e is formed after time t; &e- Since the CR dif- 
fusion coefficient k has dimensions length 2 /time, a value 
K ~ r lobe/tiobe is required to establish a CR pressure 
gradient dP c /dr over distance r; f, e in time ti Q b e . Lo- 
calized, quasi-spherical X-ray cavities must be formed 
in a time ti b e that is significantly less than the buoy- 
ancy rise time tbuoy of the cavity in the hot gas atmo- 
sphere of the galaxy group/cluster. (If ti a b e tbuoy then 
plume-like cavities would form and this is not generally 
observed.) Typically tbuoy is several times larger than 
the dynamical (freefall) time td y at the radius r where 
the cavity is formed in the group/cluster potential. Con- 
sequently, the diffusion coefficient (diffusivity) is of or- 
der k ~ 3.0 x 10 28 r 2 obe fcpc / '{Uobe/lW yrs) cm 2 s" 1 where 

~ 1 7 yrs is a typical dynamical time within a few kpc of 
the centers of large elliptical galaxies. Therefore, a diffu- 
sivity k ~ I0 30 cm 2 s _1 in the cavity region is required 
to form cavities of size 5-8 kpc. 

However, in the more acceptable calculations described 
below a somewhat smaller n ~ I0 28 cm 2 s _1 is required 
in the initial gas. The reason for this is that the cavities 
must persist for times tbuoy ~ 10 8 yrs, and for these long- 
lived cavities the relevant length scale that determines k 
is not ri b e but the thickness of the cavity wall, 5s 0.1rj o (, e - 
Consequently, cavities last for times ~ tbuoy only if the 
diffusivity in the undisturbed initial gas is n 5s 10 28 cm 2 
s _1 . This is a remarkable result since k ~ I0 28 cm 2 
s _1 is the value of the diffusivity most favored in studies 
of CR diffusion in the Milky Way disk and halo (Jones, 
1979; Strong & Moskalenko 1998; Snodin, et al. 2006). 
This noteworthy and encouraging consistency suggests 
that X-ray cavities are indeed formed by cosmic rays. At 

CR particle velocities (v) = 1.4 x 10 9 £^ cm s , the 
diffusivity k « (v)X w 10 28 cm 2 s _1 suggests a mean free 

1 /2 

path, A « 2/Ej^ ev pc, that is much less than astronomi- 



cal scales of interest, although this simple kinetic model 
for the diffusion is almost certainly inadequate. In the 
exploratory calculations described here we begin by as- 
suming a uniform diffusivity k throughout the cavity and 
gas. However, the diffusion of cosmic rays is expected to 
depend on the density of frozen-in magnetic scatterers 
which scales with the gas density. For this reason we 
also consider solutions in which k varies inversely as a 
power of the gas density, n cx p~ q . 

In solving Equation (9) it is necessary to solve the time- 
dependent diffusion equation 

at each time step over the entire grid (cavity plus gas). 
For a sufficiently large computational region we can 
safely assume that the CR flux r 2 nde c /dr vanishes at 
the outer boundary and apply reflecting boundary con- 
ditions there. However, at the center of the flow we must 
include an appropriate source of CRs. We describe the 
CR flux in the innermost computational zone using the 
steady-state point-source solution of the equation above 



where E c is the CR luminosity (ergs s x ) at the origin. 
Therefore, in solving Equation 10 as r — ► 0, we assume 
that the flux r 2 nde c /dr converges to — P c /47r, as in the 
steady state point source solution. This approximation 
should be accurate close to the CR source. 

If an X-ray cavity of radius n t, e is formed in time £; 6 e > 
then the approximate energy required to form and fill the 
cavity with CRs is 

Eiobe = C —PoVi ob e where V /o6e = -nrf obe , (12) 

7c — 1 3 

the CR luminosity is 

E c = El be/tl be, (13) 

and the characteristic cavity diffusivity is 

«o = rf obe /U be and k = c k k (14) 

where Pq is the pressure of the hot gas initially surround- 
ing the CR source and c K is a dimensionless coefficient. 
The energy E lobe includes the work done in displacing 
the gas PoViobe and the CR energy within the cavity at 
pressure equilibrium, [p/ c — 1) _1 PoViobe, which is less ac- 
curate during the period of cavity formation t 5; i; & e - 

Based on these dimensional considerations, we assume 
in the following calculations that the CR luminosity E c 
is constant for t < t/ (, e , and zero thereafter. To generate 
various types of solutions, we vary c K but keep tiobe — 
10 7 yrs and r; f, e = 6 kpc fixed in all calculations. To 
avoid an unnecessary excess of parameters, we assume 
that all cavities are produced in a uniform gas of initial 
density n e Q = 0.03 cm~ 3 and temperature kT = 2 keV 
(2.32 x 10 7 K). We assume throughout relativistic cosmic 
rays with 7 C = 4/3. With these assumptions the sound 
speed in the undisturbed gas is c s0 = 726 km s _1 , the 
initial gas pressure is Po = 1-84 x 10~ 10 dyne cm~ 2 , the 
CR luminosity is E c = 6.18 x 10 43 erg s _1 , the lobe 

energy is Ei obe = 1.94 x 10 58 ergs, and k = 1.08 x 10 30 

9 — 1 
cm' s \ 
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4. CAVITY WITH CONSTANT k 

We begin with a discussion of a solution of Equa- 
tions 5-9 in which the diffusion coefficient k = kq — 
1.08 x 10 30 cm 2 s^ 1 is constant everywhere, in the cav- 
ity and throughout the gas i.e. c K — 1. We assume 
that there is enough residual plasma in the cavity to 
sustain magnetic fields and provide inertial resistance 
to the diffusing CRs. The presence of magnetic fields 
within observed cavities is obvious from the enhanced 
radio-synchrotron emission observed there. As described 
above, the CRs are assumed to turn on at t — 0, main- 
tain a constant luminosity E Cl and turn off abruptly at 
U obe = 10 7 yrs. 

Figure 1 shows the flow evolution at four representative 
times t < 1 x 10 7 yrs. The large cosmic ray pressure gra- 
dient within about 2 kpc dominates the initial formation 
of the cavity. The leading edge of the expanding gas is 
a region of compression where the gas velocity, pressure, 
density and temperature all rise to a small maximum. 
A weak shock occurs about midway in this compression 
zone. Interior to this is a rarefaction region of approxi- 
mately uniform velocity in which the gas pressure, den- 
sity and temperature all decline adiabatically. The gas 
temperature profiles show tiny peaks with subsequent 
cooling. The flow velocities at all times shown are sub- 
sonic, u < c s o — 725 km s _1 . The small (dynamically 
insignificant) oscillations in the velocity suggest that our 
modified artificial viscosity is not entirely damping out 
the kinetic energy entering the (weak) shock within the 
leading edge of the gas flow. At time t — ti b e = 10 7 yrs, 
the cosmic rays continue to stream through the gas, pro- 
viding a pressure gradient sufficient to sustain the cavity 
in approximate hydrostatic equilibrium. 

However, just after time ti ohe = 10 7 yrs when the CR 
luminosity E c at the origin turns off, the CRs continue 
to diffuse rapidly through the gas. By time t = 3x 10 7 
years shown in Figure 2 the cosmic ray pressure has be- 
come much less than the gas pressure, allowing hot gas to 
flow back toward the center filling the cavity with nearly 
stationary gas, u « 0. The entropy near the center is 
seen to be slightly higher (lower density, higher temper- 
ature), indicating that this gas experienced a stronger 
shock when the CRs first turned on. 

The upper panel of Figure 3 shows the evolution of 
the cavity radius, r cav (t), defined (arbitrarily) as that 
radius where the gas density drops to half its original 
value. The radius rises to a maximum r cav w 3 kpc 
(about half r; (, e = 6 kpc) at t = ti f, e = 10 7 yrs, then 
drops catastrophically. Such a cavity cannot correspond 
to those observed which persist during their long buoy- 
ancy phase after the CR source has turned off. We con- 
clude that the diffusion coefficient in the gas must be 
less than kq = 10 30 cm 2 s _1 so that the cosmic rays can 
be confined within the cavity, sustaining it for a longer 
time. The lower panel in Figure 3 shows the total CR 
energy that has entered the flow E c t (solid line), the total 
CR energy in the computational grid E Citot (t) (dashed- 
dotted line) and the CR energy E CiCav (t) in the cavity 
r < r cav (dashed line). 

5. CAVITIES WITH DENSITY-DEPENDENT k 

To prolong the observable lifetime of cavities beyond 
that of the previous calculation, we consider flows in 



which the cosmic ray diffusion is slower in the gas, but re- 
tains the value no = 10 30 cm 2 s _1 in the cavity. First we 
modify the diffusion coefficient so that k — c k k where 



9(n e ) 



c K = 1 



if n e > n ed 
if n e < n ed 



(15) 



where 



g(n e ) = {n e /n ed y 



with n e d = 0.0063 cm~ 3 and q = 1.48. In the distant 
undisturbed gas c K = g{n ea ) = 0.10 is about ten times 
lower than before. However, during the initial stages 
of the calculation when cosmic rays first encounter the 
gas with a low diffusion coefficient near the origin, the 
numerical solutions can become unreliable. To avoid this 
difficulty, we also require that k have the value no near 
the origin regardless of the gas density there. To retain 
c K « 1 near the origin we adopt a final k = c K re where 



r (c K + 1) + \{c K - l)tanh 



(r- 



(16) 



and r m = Ar m = 0.3 kpc. Once the cavity radius ex- 
ceeds Ar m , which occurs very early in the calculations, 
k = kq throughout the rest of the cavity and, as the gas 
density increases, k decreases as Kog(n e ), approaching 
Kog(n e o) in the distant undisturbed gas. 

Figures 4, 5, and 6 show the evolution of this cav- 
ity. At early times t < 10 6 yrs the stronger coupling of 
CRs with the gas results in mildly supersonic velocities 
and the maximum amplitudes of the flow variables are 
somewhat greater than before. However, at later times 
t > Uobe = 10 7 yrs (Figure 5) when the central CR source 
is no longer active, the cavity again closes rapidly, nearly 
returning to its initial state by ~ 3 x 10 7 yrs. The evo- 
lution of the cavity radius shown in Figure 6 shows that 
the cavity produced at 10 7 yrs is nearly twice as large as 
that in Figure 3. However, as before, the cavity radius 
drops immediately after time ti a b e , suggesting that our 
reduced n is still too large. 

To further confine the cosmic rays, we show a final 
model similar to the previous one but with n e d = 0.0032 
cm~ 3 and q = 2.05 so g(n e o) = (n e o /n ec i)~ g = 0.01 and 
k ~ 10 28 cm 2 s _1 in the distant, undisturbed gas. Fig- 
ures 7, 8 and 9 show the corresponding cavity evolution. 
With this further reduction of CR diffusion in the gas, 
the solution at times t < M b e is similar to the previ- 
ous calculation, but with even larger amplitudes in the 
flow variables. However, Figures 8 and 9 show that the 
cavity is now able to sustain a significant size after the 
CR source has been turned off at t > ti b e — 10 7 yrs. 
The flow profiles in Figure 8 at times t = 3 x 10 7 and 
7 x 10 7 yrs illustrate this strong sustained confinement of 
CRs within r cav . The highly subsonic flow at these late 
times is nearly in hydrostatic equilibrium, but some low- 
density gas has flowed back into the cavity. The upper 
panel of Figure 9 shows that the cavity radius reaches 
a maximum of ~ 8 kpc and by 3 x 10 7 years still has 
a radius ~ 5 kpc expected from the original choice of 
parameters. 

At times t < ti b e = 10 7 yrs the CR energy in the 
cavity E c , cav and in the flow E c , to t (lower panel of Figure 
9) both slightly exceed the total energy supplied Et. The 
total CR energy is not a strictly conserved quantity and 
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since cavitation requires a net compression in the gas, 
cosmic ray energy can be created by the e c oc p lc term 
in Equation (9). However, at later times t > ti b e = 10 7 
yrs, when the cavity radius decreases, E c ^ t ~ E c t to a 
good approximation. 

The evolution of X-ray-radio cavities can also be seen 
in diffuse gamma ray emission either as inverse Comp- 
ton (IC) emission from collisions of CR electrons with the 
cosmic microwave background or as the 7r° decay contin- 
uum resulting from CR proton collisions with hot gas 
protons. It is of interest to explore how the bolometric 
gamma ray emission due to ir° decays varies with the 
evolution of X-ray cavities, assuming that protons are 
an important fraction of the total CR energy density. 
If the spectral shape of the energy distribution of CR 
protons remains approximately constant with space and 
time, the number density of CR protons is proportional 
to e c and the local emissivity from ir° decays is propor- 
tional to e c n p where n p — 0.83n e is the proton density in 
the hot gas. Using the cavity evolution shown in Figures 
7 and 8, Figure 10 shows the result of integrating e c n p 
along the line of sight, S 7 = J e c n p ds, which should be 
approximately proportional to the radial surface bright- 
ness profile of the bolometric 7r° decay continuum. 

The dramatic evolution of the pion continuum shown 
in Figure 10 suggests that it may be possible to identify 
young cavities by a strong limb brightening in the gamma 
ray surface brightness, provided k ~ 10 28 cm 2 s _1 . How- 
ever, if k < 10 28 cm 2 s _1 the pion continuum will be 
strongly peaked in the walls of cavities of all ages. In gen- 
eral the shape of the CR diffusion front in gamma rays 
provides important information about the magnitude of 
k. Since the pion continuum due to proton CRs peaks at 
70 Mev, it may be possible to detect the pion decay emis- 
sion even in the presence of an appreciable IC continuum 
due to CR electrons. Moreover, the spectral difference in 
the two gamma ray continuua could in principle allow an 
estimation of the relative numbers of CR electrons and 
protons. Finally, the level of the pion continuum at the 
center of more mature cavities can provide information 
about the density of thermal gas within the cavities and 
therefore about the rate that cavity buoyancy is trans- 
porting mass outward in the hot atmosphere. 

6. FINAL REMARKS AND CONCLUSIONS 

In these simple one-dimensional calculations we have 
shown, perhaps for the first time, how radio lobes and 
X-ray cavities are formed in galaxies and clusters. We 
have not explored a vast range of parameters but have 
concentrated instead on the critical influence of the CR 
diffusion coefficient k in cavity evolution. Of particular 
interest is how the size and duration of the cavities de- 
pend on the rate that cosmic rays diffuse into the gas. 
For the CR parameters adopted here, when the diffusion 
coefficient k is > 10 29 cm 2 s _1 , the cavities tend to be 
smaller than those observed. With these large values of 
k the cosmic rays diffuse readily through the cavity walls 
into the hot gas, rapidly establishing a quasi-steady cav- 
ity profile. But these cavities also disappear rapidly fol- 
lowing the termination of the cosmic ray source. This is 
unacceptable because cavities are observed during a long 
buoyancy time % uoy as they move radially away from 
their formation sites near the center of the group/cluster 



gas. 

However, when k SS 10 28 cm 2 s _1 in the gas, with other 
parameters unchanged, the cavities are larger, similar to 
those observed. Furthermore, such cavities persist ~ 10 
times longer than the initial formation time, an essential 
attribute of observed cavities. It is noteworthy that k ~ 
10 28 cm 2 s _1 is also the preferred diffusivity in successful 
models of CR diffusion in the Milky Way (Jones, 1979; 
Strong & Moskalenko 1998; Snodin, et al. 2006). 

Cavity formation by diffusing cosmic rays is particu- 
larly gentle since it accelerates the ambient gas over an 
extended region of thickness ~ (nt) 1 / 2 , not at a surface 
as occurs when cosmic rays are approximated with a hot 
gas. As a result the shocks produced by CRs in the 
surrounding gas are mild and rather insignificant in our 
calculations. To simulate the varying density of mag- 
netic CR scatterers that scale with the gas density, we 
consider solutions in which the diffusion coefficient varies 
inversely with some power of the gas density. Finally, the 
long-lived cavities we describe here are formed in a hot 
gas with uniform pressure and are therefore expected to 
be somewhat smaller than observed cavities which ex- 
pand as they experience lower ambient pressures during 
their buoyant outward motion in the hot gas atmosphere 
of a galaxy cluster. 

Cosmic rays are confined by partial reflection from the 
cavity walls. This point is of considerable interest since 
radio observations of X-ray cavities also indicate that 
synchrotron emitting electrons are confined within X- 
ray cavities (e.g. for Hydra A: McNamara et al. 2000; 
David et al. 2001; Allen et al 2001). Furthermore, pro- 
ton CRs with energies exceeding the tt° rest mass (135 
Mev), if trapped in cavities, will experience fewer colli- 
sions in the low-density cavity plasma, resulting in lower 
global gamma ray luminosities from 7r° — > 27 decays, 
as indicated by EGRET observations (e.g. Pfrommcr 
& Ensslin 2004), than if the CRs were uniformly dis- 
tributed throughout the hot gas as often assumed. Fig- 
ure 10 shows that 7r° production is enhanced in the dense 
gaseous walls of young cavities and in the cores of old 
cavities. Observations of the 7r° — ► 7 decay continuum 
are expected to provide important information about the 
cavity age, the CR diffusivity, the relative numbers of 
CR electrons and protons and the density and mass of 
thermal gas within the cavity. At 100 Mev GLAST is 
expected to have a spatial resolution of 30" which may 
be sufficient to resolve gamma ray profiles in the pion 
continuum in nearby groups/clusters. Inverse Compton 
emission from radio lobes is also expected to be sensitive 
to K. 

Finally, we note that during the formation of young 
cavities the cosmic ray pressure P c inside the cavity can 
exceed the external gas pressure by factors of 3-30 (e.g. 
Fig. 7) while studies of the radio emission from the cav- 
ity have traditionally assumed equal pressures within the 
cavity and the ambient gas (e.g. Dunn & Fabian 2004). 

Studies of the evolution of hot gas in elliptical galaxies 
at UC Santa Cruz are supported by NASA grants NAG 
5-8409 & ATP02-0122-0079 and NSF grant AST-0098351 
for which we are very grateful. Continued support would 
also be appreciated. 
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Fig. 1. — Flow variables for cavity evolution with constant cosmic ray diffusion coefficient k = kq = 1.08 X 10 30 cm 2 s~ 1 . The evolution 
is shown at four times 0.080, 0.183, 0.611 and 1.000 X 10 7 yrs and all profiles evolve toward larger radii as the cavity forms. In the top 
panel the dashed lines show the cosmic ray pressure and the solid lines show the gas pressure. 
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Fig. 2. — Flow variables for cavity evolution with constant cosmic ray diffusion coefficient k = kq = 1.08 X 10 30 cm 2 s 1 . The evolution 
is shown at time t = 3.0 X 10 7 yrs. In the top panel the dashed line shows the cosmic ray pressure and the solid line shows the gas pressure. 
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Fig. 3. — Evolution of cavity with constant cosmic ray diffusion coefficient. Top panel: Evolution of the radius of the cavity r cav (t). 
Lower panel: evolution of the total CR energy that has entered the flow Et (solid line), the total cosmic ray energy in the computational 
grid E c ,tot (dashcd-dottcd line) and the CR energy E c ,cav in the cavity r < r ca v (dashed line). The variation of E c ,tot closely follows E c t. 
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Fig. 4. — Flow variables for cavity evolution with density-dependent cosmic ray diffusion. The diffusion coefficient in the cavity 
k = ko = 1-08 X 10 30 cm 2 s _1 is constant but in the gas the diffusion cocfficcnt increases as k oc (ne/0.00631)" 1 ' 477 , reaching k = O.Iko 
in the distant gas where n e = n e Q. The evolution is shown at four times 0.169, 0.368, 0.668 and 1.000 X 10 7 yrs. All profiles evolve toward 
larger radii as the cavity forms. In the top panel the dashed lines show the cosmic ray pressure and the solid lines show the gas pressure. 
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Fig. 5. — Flow variables for cavity evolution with density-dependent cosmic ray diffusion. The diffusion coefficient in the cavity 
k — kq — f .08 X 10 30 cm 2 s _1 is constant but in the gas the diffusion cocfficcnt increases as k oc (ne/0.00631)" 1 ' 477 , reaching re = O.lreo 
in the distant gas where n e = n e Q. The flow profiles are shown at two times 3.00 X 10 7 yrs (light lines) and 7.00 X 10 7 yrs (heavy lines), 
fn the top panel the dashed lines show the cosmic ray pressure and the solid lines show the gas pressure. 
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Fig. 6. — Evolution of cavity with density-dependent cosmic ray diffusion coefficient, reaching k = O.Iko in the distant undisturbed gas. 
Top panel: Evolution of the radius of the cavity r cav (t). Lower panel: evolution of the total CR energy that has entered the flow Et (solid 
line), the total cosmic ray energy in the computational grid E c>to t (dashed-dotted line) and the CR energy E c ,cav in the cavity r <C 'f'cav 
(dashed line). 
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Fig. 7. — Flow variables for cavity evolution with density-dependent cosmic ray diffusion. The diffusion coefficient in the cavity 
k = Ko = 1 .08 X fO 30 cm 2 s _1 is constant but in the gas the diffusion coefficient increases as k oc (ne/0.00316)~ 2 047 , reaching k = O.OIko 
in the distant gas where n e = n e Q. The evolution is shown at four times 0.f56, 0.375, 0.600 and f .000 X 10 7 yrs. All profiles evolve toward 
larger radii as the cavity forms. In the top panel the dashed lines show the cosmic ray pressure and the solid lines show the gas pressure. 
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Fig. 8. — Flow variables for cavity evolution with density-dependent cosmic ray diffusion. The diffusion coefficient in the cavity 
k = KO = 1-08 X 10 30 cm 2 s _1 is constant but in the gas the diffusion coefficient increases as k oc (n e /0.00316)- 2 047 , reaching k = 0.01k 
in the distant gas where n e = n e o. The flow profiles are shown at two times 3.00 X 10 7 yrs (light lines) and 7.00 X fO 7 yrs (heavy lines). 
In the top panel the dashed lines show the cosmic ray pressure and the solid lines show the gas pressure. 
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Fig. 9. — Evolution of cavity with density-dependent cosmic ray diffusion coefficient, reaching k = O.OIko hi the distant undisturbed 
gas. Top panel: Evolution of the radius of the cavity r ca v(t). Lower panel: evolution of the total CR energy that has entered the flow 
Et (solid line), the total cosmic ray energy in the computational grid £ c ,tot (dashed-dotted line) and the CR energy E CyCav in the cavity 
r < r cav (dashed line). The variation of E c ,tot(t) closely follows E c t at t < 10 7 years, and closely follows E c t thereafter. 
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Fig. 10. — Approximate variation of the bolometric pion gamma ray continuum for the cavity shown in Figures 7 and 8. The gamma 
ray surface brightness S 7 is shown in arbitrary units at four times: 1.0 X 10 7 yrs (dotted line), 1.9 X 10 7 yrs (short dashed line), 1.0 X 10 7 
yrs (long dashed line), and 7.0 X 10 7 yrs (solid line). 



